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Abstract — Coherent feedback control considers purely quan- 
tum controllers in order to overcome disadvantages such as 
the acquisition of suitable quantum information, quantum 
error correction, etc. These approaches lack a systematic 
characterization of quantum realizability. Recently, a condition 
characterizing when a system described as a linear stochastic 
differential equation is quantum was developed. Such condition 
was named physical realizability, and it was developed for linear 
quantum systems satisfying the quantum harmonic oscillator 
canonical commutation relations. In this context, open two-level 
quantum systems escape the realm of the current known condi- 
tion. When compared to linear quantum system, the challenges 
in obtaining such condition for such systems radicate in that 
the evolution equation is now a bilinear quantum stochastic 
differential equation and that the commutation relations for 
such systems are dependent on the system variables. The goal 
of this paper is to provide a necessary and sufficient condition 
for the preservation of the Pauli commutation relations, as well 
as to make explicit the relationship between this condition and 
physical realizability. 

I. Introduction 

In the last twenty years, the use of quantum feedback 
control systems have became critical for the development 
of quantum and nano technologies [1], [4], [5]. However, 
the majority of approaches consider a classical controller 
in the feedback loop. In this context, coherent feedback 
control considers purely quantum controllers in order to 
overcome disadvantages such as the acquisition of suitable 
quantum information, quantum error correction, etc [3], 
[9], [14]. Unfortunately, these approaches lack a systematic 
characterization of quantum realizability. In [8], a condition 
characterizing when a system described as a linear stochastic 
differential equation is quantum was developed. Such condi- 
tion was named physical realizability, and it was developed 
specifically for linear systems satisfying the quantum har- 
monic oscillator canonical commutation relations. The class 
of systems for which this condition is known to be satisfied is 
still too limited for applications. In this paper, the focus is on 
systems describing the dynamics of open two-level quantum 
systems. Compared to a linear quantum system, the problem 
is more complicated and requires extra machinery for two 
basic reasons. The first is that the system being analyzed 
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is bilinear, and the second is that the commutation relations 
that the system has to obey are now dependent on the system 
variables, which was not the case for linear quantum systems 
related to the quantum harmonic oscillator [8], [11]. In [6], 
a characterization of the physical realizability for open two- 
level quantum systems was provided. However, it is not clear 
whether or not such condition imply the preservation in time 
of the commutation relations for the system variables of 
the bilinear quantum stochastic differential equation (QSDE) 
describing the system. Thus, the main contribution of this 
paper, given in Section [V] is to provide a necessary and 
sufficient condition for the preservation of Pauli commutation 
relations, as well as to make explicit the relationship between 
this condition and physical realizability. Furthermore, in Sec- 
tion llVI the physical realizability condition of open two-level 
quantum systems is reformulated in terms of the quadrature 
of the interacting Boson field, which yields a more natural 
self-adjoint (all component matrices are real) representation 
of the system and the physical realizability condition. 

The paper is organized as follows. Section [TT] presents the 
basic preliminaries on open quantum systems. In Section [TTTJ 
the necessary algebraic machinery to study open two-level 
quantum systems is given. This is followed by Section [IV] 
in which the definition of physical realizability is provided 
as well as a condition for a bilinear QSDE to be physically 
realizable. In Section [VJ it is shown that a physically realiz- 
able system preserves the commutation relations established 
for spin operators. Finally, Section [VI] gives the conclusions. 

II. Open Two-Level Quantum Systems 

Systems governed by the laws of quantum mechanics that 
interact with an external environment (e.g., electromagnetic 
field) are known as open quantum systems. In order to study 
such systems, one has to give a quantum description of both 
the system and the interacting environment. The quantum 
mechanical behavior of the system is based on the notions 
of observables and states. Observables represent physical 
quantities that can be measured, as self-adjoint operators 
on a complex separable Hilbert space S), while states give 
the current status of the system, as elements of f), allowing 
the computation of expected values of observables. Here 
open quantum systems are treated in the context of quantum 
stochastic processes (see [2], [13] for more information). For 
this purpose, observables may be thought as quantum random 
variables that do not in general commute. A measure of the 
non commutativity between observables is usually given by 
the commutator between operators. The commutator of two 
scalar operators x and y in S) is an antisymmetric bilinear 



operation defined as 



, y] = xy - yx. 



Also, for an n-dimensional vector of operators x and an Tri- 
dimensional vector of operators y, the commutator of x and 

y is 

[x,y] = xy T - (yx T ) T . 

In particular, the commutator of x and its adjoint x^ is the 
n x n matrix of operators 



[x, x'] = xx' — (x^x 1 Y , 



where 



„# A 



V < J 

and * denotes the operator adjoint. In the case of complex 
vectors (matrices) * denotes the complex conjugate while t 
denotes the conjugate transpose. The non-commutativity of 
observables is a fundamental difference between quantum 
systems and classical systems in which the former must sat- 
isfy certain commutation relations originating from Heisen- 
berg uncertainty principle. The environment consists of a 
collection of oscillator systems each with annihilation field 
operator w(t) and creation field operator w*(t) used for the 
annihilation and creation of quanta at point t, and commonly 
known as the boson quantum field (with parameter t). Here it 
is assumed that t is a real time parameter. The field operators 
w(t) and w*(t) satisfy commutation relations as well. That 
is, 

[w(t),w*(t')] =S(t-t'), 

for all t, t' g K, where 6(t) denotes the Dirac delta. Its 
mathematical description is given in terms of a Hilbert space 
called a Fock space. When the boson quantum field is in 
the vacuum state, i.e., no physical particles are present, it 
then represents a natural quantum extension of white noise, 
and may be described using the quantum Ito calculus [2], 
[13]. This amounts to have three interacting signals (inputs) 
in the evolution of the system: the annihilation processes 
W(t), the creation process W^{t), and the counting process 
A(t). The evolution of an open quantum system (i.e., the 
system together with the environment) is unitary. That is, if 
tp is a state then ip(t) = U(t)ip, where U(t) is unitary for 
all t, and is the solution of 

dU(t) = US - I)dA(t) + LdW\t) - L^SdW(t) 
-\^{l)L + %H) dt^j U(t), 

with initial condition U(0) = I, I denoting the identity 
operator and i being the imaginary unit. Here, *H is a fixed 
self-adjoint operator representing the Hamiltonian of the 
system, and L and S are operators determining the coupling 
of the system to the field, with S unitary. The evolution of 



(1) 



ip is equivalent to the evolution of the observable X given 
by 

X{t) = U\t){X ®I)U{t), 

whose evolution is referred as the Heisenberg picture while 
the one for ip is known as the Schrddinger picture. This 
paper exclusively takes the point of view of the Heisenberg 
picture. Quantum stochastic calculus allows then to express 
the Heisenberg picture evolution of X as 

dX = (S^XS - X) dA + C(X) dt + [X, L] dW^ 
+ [L\X]SdW, 

where C(X) is the Lindblad operator defined as 

C(X) = -i[X,H] + \ {i} [X, L] + [L\X\L) . (2) 
The output field is given by 

Y(t) = U(tfW(t)U(t), 

which amount to 

dY = Ldt + SdW. 

In summary, one can say from the discussion above that 
the dynamics of an open quantum systems is uniquely 
determined by the triple of operators (S, Hereafter, 
the operator S is assumed to be the identity operator (S — I). 

The main focus of this paper is on the dynamics of 
open two-level quantum systems interacting with one boson 
quantum field. The Hilbert space for this system is Sj = C 2 , 
the two dimensional complex vector space. The vector of 
system variables is 



Xi \ 












X2 






%3 ) 







where <7i, 02 and 03 are spin operators. Given that these 
operators are self-adjoint, the vector of operators x satisfies 
x = x#. In particular, a self-adjoint operator a in is 
spanned by the Pauli matrices [12], i.e., 



1 3 



where ccq = Tr(er), on = Tr(aai), and 



cto 



0-2 
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03 
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denote the Pauli matrices. Thus, ao and (ai,a2,«3) T G C 3 
determine uniquely the operator a with respect to a given 
basis in C 2 . The initial value of the system variables can 
be set to a;(0) = (ci, 02,173). The product of spin operators 
satisfy 



OiOj = 5ij 



k 



for k G {1,2,3}. It is then clear that the commutation 
relations for Pauli matrices are 



where 



[o~i, <Tj] = tijk&k, 



(3) 



where 5ij is the Kronecker delta and eyjt denotes the Levi- 
Civita tensor defined as 

!+l, {i,j, k} is an even permutation of {1, 2, 3}, 
— 1, {i,j, k} is an odd permutation of {1, 2, 3}, 
0, otherwise. 

Due to the fact that the Pauli matrices form a complete 
orthogonal set, any Hamiltonian and coupling operators of 
polynomial type are representable as linear functions of 
x. Therefore, assuming linearity captures a large class of 
Hamiltonian and coupling operators without much loss of 
generality, i.e., 

W = ax and L = Ax, 

where a T G R 3 and A T G C 3 . As mentioned before, the 
coupling operator specifies how the interacting field acts on 
x. In general, the dimensionality of the coupling matrix A 
depends proportionally on the number of interacting fields. 

It is customary to express QSDEs in terms of its interaction 
with quadrature fields. The quadrature fields are given by the 
transformation 



Wi 
W 2 



1 1 
— i i 



W 



(4) 



where the operators W\ and W 2 are now self-adjoint. In [7], 
the Ito table for W and is 



dW 



dW dW^ 



1 



dWi , ^00 
which in terms of the quadrature fields is 



dt. 



Observe that, in general, the evolution of x (standard form) 
falls into a class of bilinear QSDEs expressed as 

dx = F dt + Fx dt + Gix dW x + G 2 x dW 2 , (5) 

where F G M. 3 and F, Gi,G 2 G M 3x3 . The fact that all ma- 
trices in (O are real is due to the quadrature transformation 
©. The output field is 

dY = Hx dt+^ {dWi + idW 2 ) 

with H T G C 3 . Similarly, the quadrature form of the output 
fields can be obtained from the transformation 



Yx 
% 



Thus, 



dY_ x 
dY 2 



Hi 
H 2 



1 1 
—i i 



•dt 



Y 
Ft 



1 
1 



dWi 
dW 2 



(6) 



H 1= H + H* and H 2 = i{H* - H) 



are obviously real matrices. 

In this context, the goal of the paper can now be stated 
more specifically. Given a bilinear QSDE as in ©, under 
what condition there exist H and L such that (0 can be 
written as in (Q]i. Such condition is given in Section HV1 

III. Notation and Algebraic Relations 

In order to continue the description of open two-level 
quantum systems, some linear algebra identities are needed. 
Let f3 G C 3 and define the linear mapping O 

as 



-<3x3 



eos) 








-P2 


-Pz 





Pi 


P2 


-Pi 






(7) 



Note here that this definition allows j3 to be either a column 
or a row vector. The fact that p is either a column or a 
row vector will be clear from the context. It will also be 
convenient to rewrite 0(/3) in terms of its columns. That is, 



e(/3) = (e 1 (/3),e 2 (/3),e 3 (/3)) 



(8) 



The product of Pauli operators can be expressed in a compact 
matrix form thanks to the mapping 0. That is, 

xx T = I + i@(x). 

Similarly, the commutation relations for Pauli operators are 
written as 

[x,x T ] = 2iQ(x). 

Consider now the stacking operator vec : C mxn — > C rn " 
whose action on a matrix creates a column vector by stacking 
its columns below one another. With the help of vec, the 
matrix 0(/3) can be reorganized so that it gives 

/ ©i(/3) \ 
vec(0(/3)) = 2 (/3) =Ep, 

V ©3(/3) / 
where P is a column vector, 0i(/3) = efp, 



E 




and 




The set {— iei, —ie 2l — ie^} can be identified to be the adjoint 
representation of SU(2), which has as generators the Pauli 
matrices. It is thus that one can rewrite the matrix as 

fill £il2 Cil3 
£i21 £j22 £i23 
£i31 e i32 Ci33 



where the Levi-Civita tensor is known as the completely 
antisymmetric structure constant of SU (2). Observe also that 



(9) 



with i =/= j =/= k and 1^ G H ixA being an elementary 
matrix (i.e., matrix consisting of 1 in the position and 
everywhere else). In addition, the matrix E satisfies 



E T E = 21. (10) 

3 r- TO9X9 



If one defines the block matrix 1e = {lj'i}|j=i 
then E also satisfies 



EE 1 =1-1, 



and 



1f,E 



-E. 



(11) 



(12) 



The matrix lg can be identified as a tensor permutation 
matrix, which comes from the fact that the Levi-Civita tensor 
satisfies the contraction epsilon identity 



3 
i=l 



3jm&kn 3jri3kr< 



(13) 



The properties of 9(0) are summarized in the next lemma. 

Lemma 1: Let 0, 7 £ C 3 be column vectors. The mapping 
satisfies 

i. 9(0) 7 = -6( 7 )/3, 

m. 0(0)0 = 0, 
Hi. e;0(0) = fief - fid, 
iv. e(/3)6( 7 ) = 70 T - fi T jI, 

v. e (e(0)7) = [e(0), eft)], 

where / denotes the identity matrix, and ei is an element of 
the canonical basis of R 3 with i indicating the position of 
the nonzero element. 

Proof: To show one uses the fact that rj = 77 when 
77 G C. If the vector 0(0)7 is decomposed component-wise, 
then 



0(0) 7 



/3 T et7 

l T £\fi 
l T e T 2 fi 
7 T €/3 




-0( 7 )0 



Property ([m]) is true since e 



for all i and j, and 



3 

T ef = fimkfik 

k,l=l 



^ fil^ilkfik 
1 

fii^iikfik +22 Pie 



k,l=l 
k^l 

E 

fc<Z 



ilkpk 



— Pl^UkPk — ^2 Pl e ilkPk 

k<l k<l 

= 0, 

where the negative sign in the last summand was obtained 
because of the antisymmetry of euk- For property (|m[) , using 
©, it follows that 

e i e{fi) = (e i effi,e i e^fi,e i e^fi). 

The (p, q) component of the product e^ej is computed with 
the aid of dot as follows 



J 3 )pq 



E 

fc=l 



which means that e.;ej = 5yJ — Iji. Then, 

e,9(0) = (<5 a - Piei,6 i2 p - fte 2 , fe0 - fte 3 ) 
= 0(^1,^2,^3) - fid 
= fief - fid, 



which is indeed (|m[) . Property I|mj]) is obtained by applying 
(jm|). That is, 

T ei \ 

0(0)0(7) = - ( p T e 2 6(7) 
T g 3 J 
/3 T e 1 e( 7 ) 
T e 2 9( 7 ) 
/3 T e 3 9( 7 ) 

T (7^-71/) 
= -( T (7eT- 72 /) 
fi T (jeT- l3 I) 

= 7 T - V, 

Finally, the left-hand-side of (juj) can be written using (O as 

/ T (1 23 - 132)7 
9 (9(0) 7 ) = -9 T (l 3 i- li 3 )7 

\ T (ll2~ l2l)7 

The component of 9 (9(0)7) is then 



ejQ (0(0)7) e,- = 7 (l*-l <3 -)7 



= J l Jl7 ~ T l lJ 7 
= 7 T l y -0 - T l, j7 
= 7i0j - A'Tr 
Hence, from it follows that 

9 (9(0) 7 ) = 70 T ~ 07 T + 0V ~ 7 T 0/ 

= ( 7 T - 7 T 0/) - (07 T " V) 
= 9(0)9(7) -9(7)9(0) 
-[9(0), 9( 7 )]. 

■ 

The explicit computation of the vector fields in (Q]i is given 
by the next lemma. 



Lemma 2: The component coefficients of equations (01 
and are 

[x,U] = -2ie(a)x, (14a) 

[x,L] = -2iB(A)x, (14b) 

[x,lS] = -2iQ(A*)x, (14c) 

tf[x,L] = -2i6(A)A t - 2(AA t 7 - AU)i, (14d) 

[x, L f ]L = 2i6(A)A t + 2(AA 1 7 - A T A # )a;. (14e) 

Proof: For ( I14al ). one has by the definition of the commu- 
tator that 

[x, H] = [a;, ax] 

= x (ax) — ((ax) x T ) T 
= (xx ) a — (xx ) a 
= 2iQ(x)a T . 

Given that the components of a and x commute, the com- 
mutator [x, H] is rewritten in standard form by applying 
property (0) of Lemma Q] Thus, 

[x,H] = -2iQ(a)x. 

The procedure to compute (I14bl i and ( 1 14cb is identical to the 
one above. Hence, 

[x, L] = — 2i0(A) x and [x, L f ] = -2iO(A#) x. 

The computation of ( I14dl i is done by directly multiplying the 
scalar operator ]J and the vector operator [x, L}. Recall that 
a;' = x T since x is self-adjoint. It then follows that 

tf[x,L) = -2ix T A*Q(A)x 

/ A#xxi \ 
= -2i0(A) A*xx 2 

V A*xx 3 j 

= -2t0(A) (xx T ) T A^ 

= -2i0(A) (I-ie(x))A^ 

= -2»e(A)At - 2e(A)e(.T)At. 

Therefore, 

tf[x, L] = -2z6(A)A t - 2(AA t 7 - A* A) x. 

Finally, (|14et is computed similarly. That is, 

[x,lS]L = -2iQ(A*)xAx 

= -2iQ(A*)xx T A T 

= -2iO(A#)(I + iO(x))A T 

= -2iQ(A*)A T + 2Q(A*)Q(x)A T 

= 2i6(A)A t + 2(AA t 7 - A T A # )a;. 

■ 

From ( !14al i-( fT4eT ). one can now write equation © as the 
following bilinear QSDE 

dx = - 2ie(A)A f dt - 2@(a)x dt 

+ (-2AA 1 "/ + A^A + A T A # ) x dt (15) 
+ iQ(A* - A)xdW 1 -Q(A + A*)xdW 2 . 



Note that (0(A)A+)* = -0(A)At, which assures that 

i?e{0(A)A t } = i (e(A)At + (9(A)A t )*) = 0. 

Also, observe that A# — A is purely imaginary and A + A* 
is purely real. Therefore, all matrices in dT~5b are real. 

As mentioned in Section [II] the output fields Yi and Y 2 
depend linearly on L, U and the fields Wi and W 2 , i.e., 

f dY 1 \ _ f A + A# \ ( dW x \ 

\ dY 2 ) ~ { i(A#-A) ) xdt +[ d w 2 )• 

IV. Physical Realizability 

In an environment where the classical laws of physics 
apply, standard control techniques such as optimization or a 
Lyapunov procedures do not worry in general of the nature 
of the controller they synthesized. In other words, their 
implementation is always possible since the physics behind 
them still holds. However, if one desires to implement a 
controller that obeys the laws imposed by quantum mechan- 
ics (quantum coherent control), then such a task is not so 
easily achieved unless an explicit characterization of those 
laws is given in terms of the control system vector fields. 
This is exactly the purpose for introducing the concept of a 
physically realizable system in the next definition. 

Definition 1: System ® with output equation (0 is said 
to be physically realizable if there exist H. = ax, with a T € 



E. 3 , and L = Ax, with A T <G C 3 such that 

F a = -2i0(A)A t , (16a) 

F = -20(a) + AU + A T A # - 2AA t 7, (16b) 

Gx = (»(A# - A)) , (16c) 

G 2 = -0(A + A # ), (16d) 

H 1= A + A*, (16e) 

H 2 = i (A* - A) . (16f) 



Note by direct inspection that for a physically realizable 
system Gi T = — G; for i = 1,2. 

From a control perspective, it is necessary to characterize 
when a bilinear QSDE posses underlying Hamiltonian and 
coupling operators which allows to express the matrices 
comprising (0 and © as in Definition Q] Thus, the main 
result of the paper is given in the next theorem, which es- 
tablishes necessary and sufficient conditions for the physical 
realizability of a bilinear QSDE. 

Theorem 1: System (0 with output equation © is phys- 
ically realizable if and only if 

i. F = l(G 1 -iG 2 )(H 1 +iH 2 y, 

ii. G x = Q(H 2 ), 
Hi. G 2 = -Q(Hx), 

iv. F + F T + GxG\ T + G 2 G 2 T = 0. 

In which case, one can identify the matrix a defining the 
system Hamiltonian as 

a = lvcc(F - F T ) T E, 



and the coupling matrix can be identified to be 

A = ^(H 1 +iH 2 ). 

Proof: Assuming that (0 and (0 are physically realizable 
implies that (116al l- (116fT l are true. By comparison, conditions 
(|Ij)- (|mj) hold. By property (|ivj of Lemma [T] it follows that 

GiGi T = 0(A # - A) 2 

= AtA# - At A - A T A# + A T A 

- (A # A f - A # A T - AAt + AA T )7. (17) 

Similarly, 

G 2 G 2 T = -Q(A + A*) 2 

= (AA T + AAt + a # a t + A # A t)j 

- A T A - A T A# - At A - A*A# (18) 

Thus, GiGi T + G 2 G 2 T = 2 (2AA7 - At A - A T A#). One 
can now rewrite F in terms of a, Gi and G 2 as 

F = -29(a) - i (GiGi T + G 2 G 2 T ) . 

Similarly, F T = 20(a) (GiGi T + G 2 G 2 T ) since Gj = 
—G{. Hence, 

F + F T + G\G\ T + G 2 G 2 T = 0. 

Conversely, one needs to show that if conditions (|7j)- (fu>j) of 
Theorem Q] are satisfied, then there exist matrices a and A 
such that system (0 is physically realizable. Let 

<d(a) = ±(F-F T ). (19) 

It is trivial to check that the right-hand-side of (fT9l is 
antisymmetric with zero diagonal and hence this equation 
uniquely defines a via (0. Also, let A = \{H\ + iH 2 ). It 
follows that 

At A = - (HjHi +iH^H 2 - iH^Hi + H 2 ^H 2 ^ , 

A T A# = 1 (HjH* - iHTTH* + iH^H* + H* ' H*) 
and 

AAt = I (h-lHJ - iH x H 2 ] + iHzH^ + H 2 H 2 ^ . 

Recall from Section [II] that Hi = Hi* and G r = G* for 
i = 1, 2. It then follows that 

At A + A T A # - 2 AAt J 

= i (h\h x + h\h 2 - [h\H\ + h 2 h£) ij . 

From conditions and (|zzz|> . one obtains 

GiGi T + G 2 G 2 T = — GiGit — G2G2t 

= (h x H\ + H 2 H^j I - HlHi - H\H 2 . 

Therefore, 

GiGi T + G 2 G 2 T = 2 (2AA t J - A f A - A T A # ) . 



From (O, one obtains F = —F T — G\G\~—G 2 G 2 T . Then, 
0(a) - i (-2F T + GiGi T + G 2 G 2 T ) , 

which agrees with dl6b| ). Moreover, from ([§), ( fTUl ), ( fT9l and 
by applying the stacking operator to 0(a), a is explicitly 
obtained as vec(9(a)) = Ea T = -jvec (F - F T ). Multi- 
plying both sides by E T leaves 

a= lvec(F - F T ) T E, 
8 

which completes the proof. ■ 

V. Preservation of Canonical Commutation 
Relations 

The goal of this section is to show that the conditions 
presented in Theorem Q] are necessary and sufficient for 
preserving the Pauli commutation relations (0 by the system 
(0. To achieve this task, a property of the stacking operator 
and a lemma are needed. Let A £ R nXm , B £ R mxl and 
G G R lxr for n,m,l,r G N. It is well-known that the 
stacking operator used at the end of Section [III] satisfies 

vcc(ABG) = (G T ® A)vcc(B). 

Lemma 3: Let E = (ei,e 2 ,e 3 ) T , and A,B £ W lXn for 
neN. Then 
i. E T (A ® B)E = E T (B ® A)E. 
it. E T (I ® A)E = Tr(A)I - A T . 

Hi. E T (A ® B)E = A T B T + B T A T + Tt(A)Tt(B)I - 
Tr(B)A T - Tr(A)B T - Tr(AB)I. 

iv. EE T (A (g> B)E = {A® B)E + (B ig> A)E. 

Proof: To prove (0, a key observation is that the matrix 1 e 
is a symmetric permutation matrix such that 1e(A®B)1 e = 
(B®A). By O, 

E T (A®B)E = {-t E E) T {A® B)(-t E E) 
= E T (t E (A®B)t E )E 
= E T {B ® A)E. 

Expanding F T (7 ® A)E, identity is computed as 

3 

E T {I ® A)E = -^e k Ae k . 

fc=i 

From ( [T3] l. it then follows that 

fc— l / ij 

3 

= efe k Ae k ej 
fe=i 

3 / Cfcl A 

= ( e fcil7 e fci2 7 e ki3) A I e fc2 j J 

fc=i \efc3j7 



3 3 
^ ^ ^ ^ ^kir A rs € ks j 
k=l r,s=l 
3 3 

— ^ ^ Afs ^ ^ ^kir^ksj 
r,s=l fc=l 
3 

= S ' A rs (Si S S r j — Sij5 rs ) 

r,s—l 

3 

= (A T -Tr(A)l) t] 
Therefore, it is clear that 

E T (I ® A)75 = Tr(A)7 - v4 T . 
From (HJJ and equation ( fTTT i, identity dmt is obtained as 

S T (A® 5)75 
= E T {A®I){I ® B)E 
= E T {A® I){EE T + 1^) (7 $ B)E 
= (75 T (A ® I)E)(E T (1 ® 5)75) 

+ 75 T (A® 7)l s (7 ® 5)75 
= (75 T (A ® I)E)(E T (I ® 5)75) 
— E T (A ® I)\e{I ® B)\eE 
= (E T (A g I)E)(E T (I ® 5)75) - 75 T L4 ® 7)(B <g> 7)75 
= (Tr(A)7 - A T )(Tr(5)7 - 5 T ) - 75 T (AB ® 7)75 
= Tr(A)Tr(£)7 - Tr(A)B T - Tr(B)A T - Tr(AB)7 

Finally, identity ( fzTT > is obtained using dTTb and J121 . That is, 

££ T (>4 ® 5)75 = (7 - 1 B ) (A ® £)75 

= (A®5)75 + l s (,4® B)l E E, 
•> v ' 

(7? ® A) 

which completes the proof. ■ 

In order to be considered a quantum system, the system 
variables of (0 must preserve (0 for all times. The condition 
that (0 has to satisfy is 

d[x,x T ] - 2iQ(dx) = 0. (20) 

Note by the linearity of the map O that 

e(dx)= <3(F )dt + B(Fx)dt 

+ Q{G 1 x)dW l + G(G 2 x)dW 2 . 

A condition for system (0 to satisfy ( 120b is given in the next 
theorem. 

Theorem 2: Let [x t (0),Xj(0)} = 2iJ2 k eij k x k (0). Sys- 
tem (0 implies 

[x i (t),x j (t)]=2iJ2<HjkX k (t) (21) 



G\G 2 - G 2 G\ — 9(F ) = 
F T + F + G\Gi T + G 2 G 2 T = 0. 



(22b) 
(22c) 

T\\T 



Proof: From the fact that d[x,x ] = d(xx ) — (d(xx )) 
and in light of the quantum ltd formula (see [7, Theorem 
4.5]), the term d(xx T ) is expanded as 

d(xx T ) 

= (dx)x T + x{dx) T + (dx){dx) T 

= (F dt + Fx dt + Gix dWi + G 2 x dW 2 ) x T 

+ x (Fq dt + x T F T dt + x T G\ dWi + x T G^ dW 2 ) 
+ (7b dt + Fxdt + GixdWi + G 2 x dW 2 ) ■ 
(F T dt + x T F T dt + x T G\ dW x + x T G T 2 dW 2 ) 
= (F x T + xF^) dt + (Fxx T + xx T F T ) dt 

+ {G lX x T + xx T G{) dWi + (G 2 xx T + xx T G 2 ) dW 2 
+ G 1 xx T G\dW x dW x + G 1 xx T G^dW 1 dW 2 
+ G 2 xx T G\dW 2 dW x + G 2 xx T G^dW 2 dW 2 . 

The (i, j) component of d(xx T ) is computed as 

ef d(xx T )ej 

= (FoiXj + XiFQj) dt + (FiXXj + XiX T Fj) dt 
+ (G u xxj + x l x T G T 1] )dW 1 
+ (G 2i xxj + x l x T G 2j ) dW 2 
+ G li xx T Gj j dWidWx + Guxx T G 2j dW\ dW 2 
+ G 2i xx T G^ j dW 2 dW 1 + G 2i xx T G^dW 2 dW 2 

= (FoiXj + XiFQj) dt + (FiXXj + XiFjx) dt 

+ (Guxxj + XiGijx) dWi + (G 2 iXXj + X{G 2 jX) dW 2 
+ GuixdWijGi^xdWi) + G li (xdWi)G 2 j(xdW 2 ) 
+ G 2l {xdW 2 )G l] {xdW 1 ) + G 2] {xdW 2 )G 2l {xdW 2 ). 

One can compute the (j, i) component similarly. Thus, the 
(j, i) component of d([x,x T ]) is 

efd(xx T )ej — ej d(xx T )ei 

= {FiXXj — XjFix) dt + (xiFjX — FjXXi) dt 

+ (Guxxj — XjGux) dWi + (xiGijX — GijXXi) dW\ 
+ (G 2l xxj - XjG 2l x) dW 2 + (xiG 2j x - G 2j xxi) dW 2 
+ GuixdW^djixdWi) - djixdW^GuixdWi) 
+ G li (xdW 1 )G2j(xdW 2 ) - G lj {xdWi)G 2i (xdW 2 ) 
+ G^xdW^G^xdWx) ~ G 2j (xdW 2 )G li {xdW 1 ) 
+ G 2i {xdW 2 )G 2j {xdW 2 ) - G 2j {xdW 2 )G 2i (xdW 2 ) 

j] dt + [xi, Fjx] dt + [G u x, Xj] dWi 
+ [x i} Gxjx] dWx + [G 2l x, Xj] dW 2 + [xi,G 2 jx] dW 2 
+ [Gii(xdWi),Gij(xdWi)} + [G u {xdWi),G 2j {xdW 2 )] 
+ [G 2i (xdW 2 ) , Gxj {xdWx )] + [G 2i (xdW 2 ) , G 2j (xdW 2 )] 



for all t > if and only if 

d + Gj = G 2 + G\ = 



(22a) 



(Fik [xk,Xj] + F jk [xj, x k ]) dt 

-i 

n 

(G H k[xk>Xj] + G ljk [xi,Xk}) dWi 



k=l 



k=l 



(G 2 ik[xk,Xj] + G2jk[xi,x k ]) dW 2 



k=l 



V v ' 

2i e k i, n x m dt 



k,l=l 



2idki dt 



k,i=i 



J2 GukGw [x k dW 2 ,xidW 1 ] 



k,l=l 



-2iS k i dt 

n 

G 2lk G 2jl [x k dW 2 , xidW 2 ] . 

2i ^ ^ ^khnXm dt 



k,l=l 



The variation in time of the commutator [x, x T ] amounts to 

d([x, x T }) = 2i (GiG 2 - G 2 G\ + FQ(x) + Q{x)F T 
+GxQ(x)Gj + G 2 <d{x)Gl) dt 



2i(G 1 Q(x) + Q{x)Gl) dW x 
2i (G 2 Q(x) + 9(x)Gf ) dW 2 . 



(23) 



Replacing (l23l into ( f20l > amounts to 

2i {G 2 Gl - G X G\ - 6(F ) + FQ(x) + G(x)F T 
+G 1 Q(x)Gj + G 2 Q{x)G^ - @(Fx)) dt 
+ 2i (Gi6(x) + e(x)Gf - 0(Gix)) dW x 
+ 2i (G 2 Q(x) + Q(x)G 2 - Q{G 2 x)) dW 2 = 0. (24) 

From [13, Proposition 23.7], one can also equate the inte- 
grands in d24l i to zero. Thus, the equations to be satisfied for 
preservation of commutation relations are 

G 1 Q(x) + &{x)Gl - Q{G lX ) = (25a) 

G 2 0(x) + Q(x)G^ - <d(G 2 x) = (25b) 
G 2 G\ - G\G\ - 9(F ) + FQ(x) + S(x)F T 

+ G 1 Q(x)Gj + G 2 Q{x)Gl - Q(Fx) = 0. (25c) 



For ( 125 ai and (I25bl i. apply the operator vcc and multiply by 
E T to the left 

E T (I ® Gi)Ex + E T (Gi ® I)Ex - E T EG l x = 0. (26) 

From identities (01 and (Imb of Lemma [3] and ( [Tol l, 

(Tr(Gi)J - G t T -Gi)x = 0. (27) 

Similarly for (I25cb . one has that 

vec(G 2 Gf - GiG^ - 9(F )) 
+ (I ig> F)Fx + (F (g> 7)Fx - FFx 
+ (Gi ® G{)Ex + (G 2 <8> G 2 )F:r = 0. 

From identities dH]i and (Imb of Lemma [3] and (fTOb . 

F T vec(G 2 Gf - G X G T 2 - 9(F )) + 2 (Tr(F) - F t - F 



J 1 = 0, 



A key observation is that x(0) is represented by the linearly 
independent Pauli matrices, and that any linear combination 
aiXi(0) +a 2 x 2 (0) +03X3(0) 7^ unless 01 = a 2 = a 3 = 
for ai , a 2 , 03 G C. In addition, no linear combination of Pauli 
matrices generates the identity. So, given that x(0) ^ 0, any 
equation involving the system variables of the form Ax = b 
(A 6 C 3x3 ) implies A and b must be identically 0. Thus, 

Tr(Gi)/-Gf -Gi =0 (28a) 
Tr(G 2 )I - G 2 - G 2 = (28b) 
G X G T 2 - G 2 G\ - 9(F ) = (28c) 

, (Tr(G,:)) 2 



Tr(Gf) 



Tr(F)/-F T -F + ^((Gf 
- Tr (G,)Gf - 



-I 



7 =0. 



(28d) 



The trace of G\ and G 2 can be calculated from ( I28ab and 
(l28bl . That is, 

Tr(G 4 + Gj - Tr(G t )/) 

= Tr(Gj) + Tr(Gf ) - Tr(G ?; )Tr(7) 
= 2Tr(G,) - 3Tr(Gi) 
= -Tr(Gi) = 0, 

which leaves Tr(G 4 ) = 0. Thus, Gi = — Gj. Applying this 
result to ( I28db gives 



F T + F-Tr(F)I = j2((Gf) 2 -^fll 



(29) 



Similarly, applying the trace operator to the previous equa- 
tion gives 

2Tr(F) - 3Tr(F) = £ (tt ((Gj) 2 ) 

which amounts to 

Tr(F) = i(Tr(G?)+Tr(G2)). (30) 

Replacing (f30b into ( 1291 and since G, = — G 4 T , one has that 

F T + F + GiGi T + G 2 G 2 T = 0. 

Conversely, it is going to be shown that if d22ab - d22cb hold 
then ( |24T i holds as well. One can see from d22a| ) and ((m]) of 
Lemma [3] that 

(Tr(Gi) - Gf - G») = (E T {I ® G { )E - G t ) = 0. 

Applying $vvh of Lemma [3] and multiplying on the left by 
F and on the right by x gives 

EE T {I ®Gi)Ex - Edx 

= {I® Gi)Ex + (Gi ® 7)Fx - EGix. 



Let vec -1 denote the inverse of the stacking operator. That 
is, vec _1 (vec(yl)) = A for an arbitrary square matrix A. 
From the definition of 9, it follows that 

vec" 1 ((/ <g> Gi)Ex + (Gi <g> I)Ex - EG.x) 

= vec" 1 ((7 (g> G 4 )vec(9(x)) + (G, ® l)vec(0(x)) 

-vcc(Q(G z x))) 
= Gi@{x) + B(x)Gj - Q(G iX ). 

Thus, 

Gi@(x) + <d(x) T Gf - Q(Gix) = 0, i = l,2. (31) 

Note that (I22bb appears explicitly in the first line of d24l i. 
Next, computing the trace on both sides of d22cl i gives 

2Tr(F) = Tr (Gf) + Tr (G 2 .) (32) 

Equations fl32l and ( |22at allow one to write J22ct as 

^ + F - Tr(F)7 - £ ((Gf) 2 + <?^J 

-TrCG^Gf-^Ml/^O. 

From and ( Imt of Lemma [3] it follows that 

i 2 

F - E T (I ® F)F - - ^ E T (Gi ® G»)-E = 0. 
" i=i 

By identity dwl i of Lemma [3] and multiplying on the left by 
F and on the right by x, (122d i is equivalent to 

EFx -(I® F)Ex - (F ® I)Ex 

- (Gi ® Gi)Ex - (G 2 ® G 2 )Ex = 0. 

Applying vec -1 gives 

F9(i) + 0(x)F T - O(Fx) 

- G 1 Q(x)Gl + G 2 <d{x)Gl = 0. (33) 

Obviously if OTb and ( T33l > hold then d24l i is zero, which 
completes the proof. ■ 

Theorem 3: A physically realizable system satisfies the 
conditions of Theorem |2] 

Proof: It is enough to show that if conditions ([!]), ([mJ and 
of Theorem Q] hold then (I22at - (l22cb are satisfied as well. 
Define H = \{H\ +iH2). By conditions ([m]) and (|zzz[> of 
Theorem Q] one has that 

F Q = ^(G 1 -iG 2 )(H 1 +iH 2 )l = -2iO(H)Hl 

Now, from Lemma [T] property (jwj), 8(^0) is 

0(F O ) = -6(2i9(/f)/f t ) 

= -2i (Q{H)e(H^) - 6(i? t )e(i/)) . 



Since G; = -Gf , it follows that 

G 2 G\ - G x Gl = -2i{Q{H)@(H r ) - Q(H*)Q(H)). 

Therefore, j22bj holds. Again, from (Jm|) and (|mz|) of The- 
orem [1] a physically realizable system satisfy Tr(G;) = 
and Gj = —Gi, which imply condition d22a| ). Finally, it is 
clear that condition ( 122cb is equivalent to 

F T + F + GiGi T + G 2 G 2 T = 0, 

which concludes the proof. ■ 

VI. Conclusions 

A condition for physical realizability was given for open 
two-level quantum systems. Under this condition it was 
shown that there exist operators TL and L such that the 
bilinear QSDE (f5]l with output equation (|6]l can be written as 
in Also, it was shown that physical realizability implies 
preservation of the Pauli commutation relations for all times. 
Future work includes extending the formalism for the case 
of multi-particle spin systems. 
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